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ABSTRACT 


Large space structures are difficult to control because of the high order of their 
mathematical models. The high order mathematical model makes the use of a reduced 
order model to control the structure desirable. The Karhunen-Loeve expansion along 
with Galerkin’s method is used to generate a reduced order model. A control algorithm 
is achieved by applying linear quadratic regulator theory to the reduced order model. 

The Karhunen-Loeve basis functions or mode shapes must first be found to identify 
the reduced order model. Previous results have shown that in the limit as the structural 
damping approaches zero the Karhunen-Loeve mode shapes and natural mode shapes 
converge. Numerical techniques are applied to evaluate the structural damping required 
for convergence. Once the Karhunen-Loeve mode shapes are determined, the reduced 
order control model is applied to the full order system. The performance of various 
Karhunen-Loeve models is compared by measuring the modal energies in the controlled 


and uncontrolled modes. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. LARGE SPACE STRUCTURES 

The lightweight, flexible materials used to construct large space structures (LSS), 
like a space station, are lightly damped and when disturbed will vibrate for a consider- 
able amount of time. This prolonged vibration could jeopardize the structural integrity 
of the structure or disturb experiments on the LSS. The purpose of this thesis is to studv 
the effects of controlling this structural vibration with a reduced order model, specifically 


a Karhunen-Loeve reduced order model. 


B. PROBLEM APPROACH 

The solution of the vibration control problem requires the use of a mathematical 
model that describes the behavior of the system in time. The space structure is modeled 
as a combination of small plates of unit mass connected to form the complete structure. 
The vibrational motion of these plates can be modeled as a set of coupled damped har- 
monic oscillators. Using modal analysis the mathematical model of the structure can 
be decoupled to yield a set of uncoupled simultaneous second order differential 
equations. 

For LSS this model 1s of very high order, making control design difficult. It 1s 
therefore desirous to use a reduced order model (ROM) to control the structure. A 
control system 1s designed for the LSS using the Karhunen-Loeve (reduced order) model. 
The control system is then applied to the LSS after it has been disturbed by an impulse. 
The performance of this system and a control system based on a modal model, truncated 
by frequency, are compared. 

The LSS used in this thesis is a dual keel space station, see Figure 1 on page 2. A 
mathematical model of this space station is provided courtesy of McDonnel Douglas 
Astronautics. A computer simulation of this space station is used to determine the ef- 
fectiveness of the control system. Karhunen-Loeve (reduced order) models of increasing 
size are simulated and the modal energies are calculated and used to determine the rela- 
tive effectiveness of the control models. These results will then be compared to results 


obtained for the modal model truncated by frequency. 
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Figure 1. Representation of a Dual Keel Space Station from [Ref. 1: p. 1] 


C. ORGANIZATION 


Chapter II 1s a statement of the mathematical model of the space station. 


Chapter III discusses the Karhunen—Loeve expansion, the relationship between 
the Karhunen—Loeve mode shapes and the natural mode shapes and how a re- 
duced order model can be synthesized using the Karhunen—Loeve expansion and 


Galerkin’s method. 

Chapter IV discusses the determination of the Karhunen—Loeve modes. 

Chapter V describes the control solution and presents the control simulation re- 
sults. 


Chapter VI discusses the conclusions based on the simulation results and a re- 
commended design procedure. 


Il. THE MATHEMATICAL MODEL 


A. INTRODUCTION 

Large space structures are flexible, lightly damped structures that can vibrate for a 
considerable amount of time when disturbed by an external force. To effectively control 
such a system, a mathematical model describing the evolution of that system in time 1s 


required. 


B. THE MODAL MODEL 

The space station is modeled as a combination of small plates of unit mass con- 
nected to form the complete structure. This model can be visualized as a svstem of 
masses connected by springs and dashpots (the springs representing the stiffness factor 
and the dashpots representing the damping factor). The displacement of the masses can 
be described by a second order matrix differential equation of motion: 


M(t) + = Kx(t) + Kx(t) = F(t) (1) 
where 
e xis the generalized coordinate vector 


e Ms the diagonal system mass matrix 


Gr. 
e wo, K is the structural damping term 
e duis the damping constant 
e w,is the frequency of oscillation of the system 
e kK is the symmetric system stiffness matrix 


e F(t) is the system forcing function 


This equation represents a system of simultaneous, second order differential equations 
that are coupled by the K matrix.[Ref. 1: p. 3] 

This equation can be uncoupled and the system represented bv a set of independent 
second order differential equations. This is done through the process of modal analysis 
which is outlined in [Refs. 2, 3]. The resulting modal model consists of a set of inde- 


pendent second order differential equations: 


[#+dQyn+Qy =X’ FI (2) 


where 
°® » 1s the coordinate vector or modal amplitude vector 
e (Y= diaglw,, @,. ..- ,@,_] 
e X’=[x,x, ... x,J the transpose of the modal matrix or mode shape vector 


e F 1s the torquing force applied at a point 


Next a discrete-time equation describing the motion of the space station in terms of its 
natural modes of vibration is developed. The discrete-time state equation for the ith 


equation of motion is: 
Z(AT + 1) = O(T) Z(KT) + (1) X/ (FAT) + WD] (3) 


where 


e Z,1s a vector of the ith modal amplitude and 
the 1th modal velocity 


e , 1s the ith state transition matrix 

e J7,1s the ith input vector 

e X/ is the transpose of the ith mode shape vector 

e F 1s the control torque force vector applied at a point 
e TT 1s the sampling time 

e kis the time index 


e W 1s the disturbance input 


This equation is used for computer simulation of the space station and control 
solution.[Ref. 1: p. 4] 


Il. APPLICATION OF THE KARHUNEN-LOEVE EXPANSION TO THE 
REDUCED ORDER CONTROL OF LARGE SPACE STRUCTURES 


A. INTRODUCTION 

Large flexible structures, such as a space station, as a class of distributed parameter 
systems (DPS) require a finite dimensional model for control design. This model may 
be achieved by approximating the state of the LSS using the Karhunen-Loeve (KL) ex- 
pansion. The expansion is truncated to provide the finite dimensional approximation 
of the state for control design. The KL model that results describes the evolution of the 
approximated state of the structure. 

The natural mode shapes are normally used for modeling and control of flexible 
Structures. The relationship between the natural mode shapes and the KL mode shapes 
is described in this chapter as well as the use of Galerkin’s method. Galerkin’s method 
is used to generate a reduced order model (ROM), using both the natural mode shapes 


and the KL expansion. 


B. THE KARHUNEN-LOEVE EXPANSION 

The purpose of the KL expansion 1s stated by Stark and Woods [Ref. 4: p. 322]. 
“The idea [of the KL expansion] 1s to decompose a general second-order random proc- 
ess into an orthonormal expansion whose coefficients are uncorrelated random 
variables.” The state of a large space structure (LSS), y(x), can be modeled as random 
process since it depends on random excitations, 1.e., noise from onboard machinery and 
actuators. The second order moments of the LSS are proportional to the physical en- 
ergy and therefore exist. The LSS can therefore be approximated using the KL expan- 
sion [Ref. 5: p. 12]. This is done by projecting the random process onto an orthonormal 
basis and truncating to \ terms. The value chosen for N is a matter of “engineering 
jeeement . (Ref. 5: p. 13} 


The selection of the orthonormal basis is made by solving the KL eigenequation: 
= yess), Pz) = Ai (x) (4) 


where 
¢ (x) is referred to as an eigenfunction (or the KL mode shapes) 


e 4,1s the eigenvalue and is a measure of the excitation of the ith basis function 


e =R,,(x,z) 1s the correlation function of y(x): R,,(x,z) = ELy(x) y7(x)] 
e <e, +> isan inner product: <a(z), b(z) >= J ar(2) M(z) b(z) dz 
e M(z) is the mass density of the structure 


e is the spatial extent of the structure 


The state of the LSS 1s approximated by: 


re 
va = > Sb Ax) (5) 
j=l 
where 
® y,1S an approximation to the state of the space structure 
e () is a setuol coordinatessiouncsp — <a (a ee 


e (x) is the ith basis function or KL mode shapes 


The expansion is truncated by keeping the eigenfunctions (KL mode shapes) associated 
with the N largest eigenvalues. [Ref. 5: p. 13] 

The KL expansion vields the best approximation to the random process, 1.e., mini- 
mizes the expected value of the norm of the error, of any orthogonal expansion. The 


approximation error 1s defined as: 


e(x) =x) —yalx) = >) Gb) (6) 


For a proof of the optimality of the KU expansion see [Ref. 6:p. 11)"or | Ref. 52 pay 


C. RELATIONSHIP BETWEEN THE KARHUNEN-LOEVE MODE SHAPES AND 
THE NATURAL MODE SHAPES 

The relationship between the KL mode shapes or KL basis functions and the natural 
mode shapes is taken from Burl [Ref: 5: p. 13]. The KL mode shapes of a structure are 
related to the natural mode shapes of that structure by a linear transformation repres- 


enting a change of basis which can be written: 
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where 


r, fn) 80 me) 0 mx). | 
i @=| Q (x) QO = (x) 0 °) 
c| = Le c ¢ et | (9) 


and {7{x)} is the set of natural mode shapes. The state of a structure consists of a 


generalized position and velocity which can be expanded in terms of the natural mode 


shapes: 
y(x,1) = a4) “A” | + 400) nto | =n" (x)a(0) (10) 
j=) 
where 
a ()=Lon(t) a(2) o(0) b(2) (1)... J (11) 


are the coordinates and velocities of the natural mode shapes. The vectors, c, can be 


found by solving the equation: 


EL a(t) «(1)] Qe =A; ¢; (12) 


2 2 2 
dentin: WW () We () W56 () 


Equations 10 and 11 give the KL mode shapes in terms of the natural mode shapes. 


where 


Equation 12 is an infinite dimensional eigenvalue problem that can be solved practically 


by truncating it to the most significant terms. 


D. REDUCED ORDER MODELING USING THE KARHUNEN-LOEVE MODES 
Galerkin’s method is used to produce reduced order state equations. This can be 
done with either the KL modes or the natural modes. 


The discrete-time state equation of a distributed parameter system is: 
Wx,k) = Fy(x,k — 1) + G f{k) (1a) 


where 
® y(x,k) is the state 
e f(k) 1s the input 
e kis the time index 
e Fis an operator on the state space 
¢ Gis an operator from &” to the state space 


e For each k, v(x,k) e the state space, f(k) e 2” 


A finite dimensional approximation to this equation can be obtained using Galerkin’s 
method 


P, y(x,k) = (P,, F P,) y(x,.k — 1) + P,, G u(k) (15) 
where P, can be written in terms of a basis [y, | i= 1,2,..., 7] 
n 
P,(©) =) nlx) < nfo), +> (16) 
i=] 


The KL mode shapes (basis functions) can be used in the above equations and the KL 
model results. The natural mode shapes can be used in the above equations, which 1s 
equivalent to truncating the modal equations, producing the modal model.[Ref. 5: p. 
13) 
E. SUMMARY 

The KL expansion can be used to approximate the state of a LSS. This is done by 
projection of the state, which is modeled as a random process, onto an orthonormal 
basis function. The basis function can be found by solving the KL eigenequation. These 
KL basis functions will yield the best approximation to the state of the LSS. Then using 
Galerkin’s method a reduced order model of the LSS 1s produced. This reduced order 


model is used to generate a control which 1s applied to the entire system. 


IV. DETERMINATION OF THE KARHUNEN-LOEVE BASIS 
FUNCTIONS 


A. INTRODUCTION 

The KL basis functions (or mode shapes) can be determined bv solving the 
eigenequation 12. This is a tedious and laborious processes. It can be shown that for 
lightly damped structures the KL mode shapes can be determined from the open loop 
response by ordering the natural mode shapes in order of decreasing modal energies. 

The KL mode shapes were calculated using the KL mode program in Appendix B. 
These calculations were compared to the mode shapes selected bv observing the open 
loop response to verify that this 1s a valid method of determining the KL mode shapes. 

The value of the structure’s damping factor used in the program was increased until 
the KL mode shapes no longer were the same as those determined from the open loop 
response. This was done to determine how large the structure’s damping factor could 


be and still have the KL mode shapes converge to the natural mode shapes. 


B. NUMERICAL DETERMINATION OF THE KARHUNEN-LOEVE MODE 
SHAPES 

The KL mode shapes are found bv determining the eigenvalues and eigenvectors 
(equation 12) of the covariance matrix Elazr)«(1)]. This is an infinite dimensional 
matrix which is truncated to a finite dimensional square matrix. There are three terms 
in El.a(r) «7(r)], where a(z) is defined by equation 1]. The first of these terms is com- 


puted, for a white noise input 


ELa,(t) «(t)] = ro h(t)ot (17) 


where h(t) is the impulse response for a,(t) [Ref. 6: p. 66]. The impulse response is given 
by: 


” 


b,. 
h(t) = $$ ail SIn{Wo/ 1 — y t) (18) 


Doin) = y 


where 


¢ b,,is the modal slope in the x direction 


¢ ys the damping coefficient; y = £ iS a constant 


¢ w,,1s the natural frequency 


Performing the integration in equation 17 produces the result: 


ELat) (01) = Ki K,| 2 _ ——4____ | 19 
od 71 2a? +b] ieee a7) 
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The other terms in E[ a(t) «7(t)_] can be found in a similar way. They are: 
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Equations 19, 20, and 21 specify the eigenvalue problem equation 12. 
The eigensolution specifies the transformation from the natural mode shapes to the 
KL mode shapes. This finite dimensional eigenvalue problem is solved and the KL mode 


shapes are approximated as a linear combination of the first fiftv (flexible) natural mode 


shapes. The program in Appendix B computes the KL mode shapes by solving the 


eigenvalue problem equation 12. 


C. EMPIRICAL DETERMINATION OF THE KARHUNEN-LOEVE MODE 
SHAPES 

The numerical computation of the KL mode shapes 1s difficult. Solving the 
eigenvalue problem requires the determination of the impulse response and solving for 
all the terms in the covariance matrix EL@,(t) o,(t) J. It turns out that, in the limit as the 
damping factor of the structure approaches zero, the KL mode shapes converge to the 
natural mode shapes [Ref. 7]. 

If the modal energies for the open loop response are computed and the mode shapes 
ordered by decreasing modal energy, these mode shapes should be the same as those 


computed numerically. The output of the system, y, 1s defined: 


7 9 0 0 
7 9 1°? 9 Oy 
ei i tall (22) 
56 0) 4) eee Ws56 () X56 
0 0 0 | 
where y, € & and 
E[ith Modal Energy] = £T |l 3; 171. (23) 
The energy, given in equation 23 can be written [Ref. §& }: 
ELA IPI = | “WrydolPe (24) 
0 


where h,, is the response of y, due to an impulse applied at the disturbance input (node 
69 or 55). Equation 24 is evaluated using computer simulation. The mode shapes found 
are compared to those determined numerically for noise input at node 69 or node 55, see 


Table 1] on page 12. 


Table 1. DETERMINATION OF THE FIRST FIFTY FLEXIBLE 
KARHUNEN-LOEVE MODE SHAPES 
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Table I shows that the mode shapes are the same. Therefore, the KL mode shapes can 
be determined from the open loop response when the structural damping factor is suffi- 
ciently small. In the case of the space station simulated in this thesis a value of 0.001 
was used. The next question is. how big can the damping factor of the structure be be- 
fore the KL mode shapes fail to converge to the natural mode shapes? 
1. Required Magnitude of the Damping Factor for Convergence 
When the structural damping factor is small enough, the eigenvectors solved for 


in equation 12 approximate the natural basis 
ii : 
or 
e, =[0100...] (26) 


etc.; one element is unity and the other elements are zero. When these vectors are sub- 
stituted into equation 7 it is easy to see that the KL mode shapes are the natural mode 
Shapes. If the damping factor is increased, at some value the eigenvectors have more 
than one non-zero element and the KL mode shapes become a linear combination of 
natural mode shapes. 

The damping factor was increased successively by a factor of two from a starting 
value of 0.0005. The norm of the error (or error norm) is used as a means of measuring 
how closely the eigenvectors obtained from the KL mode program approach the ideal 


of equation 25. The error norm 1s defined as follows: 






50 
| € li, = minimum| oy — ¢; || = U-xf+ ) x (27) 
if] 
=2 
where 
ch = [xX, X)X3 ... Xj + Xs] (28) 


and for simplicity x, is assumed to be positive and the largest component of c,. Equation 


27 can be simplified to 


13 


he lk = (29) 





Since the eigenvectors are normalized, 


and substituting equation 30 into equation 29 yields 
let = Y2 (1 — x) (31) 


a. Results 
The norm of the error, ¢, and the natural mode shape sequence obtained 
from the KL mode program were used to determine how large the damping factor could 
be and still have the KL mode shapes and natural mode shapes converge. The error 
norm is presented graphicallv for each KL mode; all cases are for disturbance torques 


due to actuator noise. The cases presented are: 


¢ Figure 2 on page 16; for a damping factor of 0.0005 all but two KL modes have 
an error norm below 0.2. The natural mode sequence 1s as in Table | on page 12. 


e Figure 3 on page 17; for a damping factor of 0.001 nearly 74 percent of the KL 
modes have an error norm of 0.2 or below, the natural mode sequence 1s presented 
in Table 1 on page 12. The natural mode shapes are a good approximation of the 
KL mode shapes. 


e Figure 4 on page 18; for a damping factor of 0.002 onlv 46 percent of the KL 
modes have an error norm at or below 0.2, and the natural mode sequence no 
longer conforms to that in Table | on page 12, 1.e., the KL modes and natural 
modes are not converging. 


e Figure 5 on page 19; for a damping factor of 0.005 only 28 percent of the KL 
modes have an error norm of 0.2 or less and the KL modes and the natural modes 
are more divergent. 

b. Conclusion 
If the damping factor is greater than 0.001 the natural mode shapes are not 
a good approximation of the KL mode shapes. The convergence criteria 1s that the 
norm of the error, for 50 percent of the KL modes or greater, is more than 0.2. The 


value for the damping factor used in the space station simulations as noted in previous 


14 


chapters was 0.001. The assumption that the KL mode shapes and natural mode shapes 
converge is valid for this damping factor and determining the KL modes from the modal 


energies 1S a good approximation. 
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Figure 2. Norm of the error values for the KL modes, excitation at node 69 , 


damping factor = 0.0005 
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Figure 3. Norm of the error values for the KL modes, excitation at node 69 , 


damping factor = 0.001 
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Figure 4. Norm of the error values for the KL modes, excitation at node 69 , 


damping factor = 0.002 
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Figure 5. Norm of the error values for the KL modes, excitation at node 69 , 


damping factor = 0.005 


VY. CONTROL SOLUTION 


A. INTRODUCTION 

The response of the space station to disturbance and control inputs is simulated 
using the modal model discussed in Chapter II]. The Fortran program that simulates this 
model was written by Preston [Ref. I: p. 47], and is used in this thesis with minor mod- 
ification, (see Appendix A). The objective of the simulation is to determine the svstem 
response to disturbances applied at modes 69 and 55 (see Figure 6) on the structure, 
using increasing sizes of KL models (1.e., 5, 10, 20 modes). The response of the system 


is depicted graphically by displaving the energy in each mode in english units of inch— 


pounds (in—lbs). 





Figure 6. Disturbance Locations 
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B. FULL ORDER SPACE STATION MODEL 

The plant being controlled in this thesis, as stated in Chapter I, is a preliminarv 
version of the NASA dual Keel space station. NASTRAN was used to generate the first 
fiftv flexible modes (starting with mode seven). The modal model that forms the full 


order model is composed of these flexible modes: 


0) 1 0 () 
X4 — 05 es dw, eee 0) () x7 
X56 0 0 0 ] X56 
0 0 wy, — s¢ 
(2) 
band Q7 (on 
u 
Ls | 
we eas 
where 
__ | modal amplitude By 33 
“11 modal velocity G3) 
and 


®¢ w, 1s the natural frequency of the ith mode 
¢ d= 0.001 is the damping coefficient 


e ue & is the control input torques from three orthogonally oriented control mo- 
ment gyros 


e we # is a random, white noise disturbance input (three orthogonal disturbance 
torques) 


e op) € & is the modal slopes at p 
¢ p,is the location of the control moment gyros (node 69 on Figure 6 on page 20) 


e p, is the location of the disturbance input (either the control moment gyro location, 
node 69, or the alpha joint, node 55, as seen in Figure 6 on page 20). 


2\ 


This model is used as the full order model in the simulation of the space station. How- 
ever, it 1s itself a reduced order model since the “infinite” number of vibrational modes 
of the structure are truncated to fifty. 

It is assumed that perfect modal amplitude and modal velocity information is avail- 
able. This assumption of perfect sensor information 1s not realistic, but isolates the ef- 


fect of modal truncation on control algorithm synthesis. 


C. THE REDUCED ORDER MODEL 

The reduced order model of the space station is obtained by truncating modes. The 
criteria for truncating modes is either the modal frequency or the Karhunen-Loeve or- 
dering (the KL ordering found in the previous chapter to be equivalent to an ordering 


based on the energy measured during open loop excitation). 


D. PERFORMANCE FUNCTION AND OPTIMAL CONTROL 

The full-order model and a reduced order model have been established. Next, a 
mathematical expression of svstem performance is required. The performance function, 
J, stated here is developed in detail in [Ref. 1: p. 10]. J consists of the total energy in 


the modeled modes plus a control energy term: 


a; 0 0 9 
() ] @@@ 0 () x7 
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where 

¢ imax is the number of modes in the reduced order model 

e r=10-" is a weighting coefficient on the control energy term selected to yield time 

constants on controlled modes of approximately 30 seconds. 

Vibration damping is achieved by application of steady state, linear quadratic regulator 
theory [Ref. 9] to the KL (reduced order) model. The control torque vector, u(k), is the 
product of an optimal gain matrix, L, and the time varying state matrix, Z (equation 3). 
The L matrix is found by solution of the Ricatti equations to minimize the performance 
function, J, [Rett pare 
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E. SIMULATION RESULTS 

The space station was simulated and the modal energies calculated for disturbance 
inputs applied at nodes 69 and 55 as shown 1n Figure 6 on page 20. The control system 
consists of KL (reduced order) models of five, ten or twentv modes. The data is pre- 
sented graphically as energy in each of the modes. The open loop response 1s included 
for comparison except where the difference in scale precludes it. The closed loop re- 
sponse for a reduced order model generated by modal truncation 1s included for com- 
parison with the KL model. The cases presented are: 


¢ Open loop (no control) response which identifies the reduced order model to be 
msed. See Figure 7 and Figure 8. 


¢ Closed loop response for a reduced order model generated by modal truncation 
with the modes ordered by natural frequency. See Figure 9 to Figure 14. 


¢ Closed loop response for a Karhunen-Loeve (reduced order) model. The KL modes 
Were determined from the open loop responses. These results are shown in 
Figure 15 through Figure 19. 


¢ Closed loop response, with a KL model based on the open loop response, for the 
system excited by the control actuators. These results are shown in Figure 20 
through Figure 22. 


These results are discussed in detail in Chapter VI. 
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Figure 7. Open Loop Response, excitation at node 55. 
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Figure 9. Five Controlled modes selected by frequency, excitation at node 53. 
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Figure 10. Ten controlled modes selected by frequency, excitation at node 55. 


NODE 90 ENERGY PER MODE 


—4 


*10 


CUO ae-0 co.0 


LEGEND 
UNCONTROLLED 


20 CONTROLLED MODES 





ENERGY(IN—LBS) 
7.5 10.0 12.5 15.0 17.5 


0.0 


Oo 


C EEE | dk de G8 
o af qin! oly ie A Ss 


6 9 121518 212427303336 39 42 45 48 51 54 57 60 
MODE 


Figure 11. Twenty controlled modes selected by frequency,excitation at node 55. 


28 


NODE 69 ENERGY PER MODE 





1s S 
a 
NS 
© 
LEGEND 
off _ UNCONTROLLED 
s| —5 CONTROLLED MODES 
ne 
8 
| 
— 
CO © 
Po os 
ae 
Z 
Pa 1c 
: | 
a i} gd pe 


6 9 1215 18 212427303336 39 42 45 48 51 54 57 60 
MODE 


Figure 12. Five controlled modes selected by frequency, excitation at node 69. 
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Figure 13. Ten controlled modes selected by frequency, excitation at node 69. 
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Figure 14. Twenty controlled modes selected by frequency, excitation at node 69. 
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Figure 15. Five controlled modes selected by modal energy, excitation at node 55. 


a2 


NODE 355 ENERGY PER MODE 


0.030 


LEGEND 
UNCONTROLLED 
10 CONTROLLED MODES 


0.025 





eS 
CJ) 
oO 
= 
ee) 
m 
uD 
a3 
Se 
as 
rr 
=e 
[x] © 
eS 
uD 
CS 
© 
ce) 
eS 
ce) 
S : 
© jn. = nee 
6 9121518 212427303336 39 42 45 48 51 54 57 60 
MODE 


Figure 16. Ten controlled modes selected by modal energy, excitation at node 55. 
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Figure 17. Twenty controlled modes selected by modal energy, excitation at node 
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Figure 18. ‘Five and ten controlled modes selected by modal energy, excitation at 
node 69. 
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Figure 20. Five controlled modes selected by modal energy due to excitation at node 


69, actual excitation at node 535. 
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Ten controlled modes selected by modal energy due to excitation at node 


69, actual excitation at node 55. 
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Figure 22. Twenty controlled modes selected by modal energy due to excitation at 


node 69, actual excitation at node 55. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Results supported theoretical claim that the KL mode shapes converge to the na- 
tural mode shapes. For the space station example, a damping factor of 0.001 was re- 
quired to make the approximation based on convergence good. 

Preston discovered, see [Ref. 1: p. 46], that truncating the modes based on natural 
frequency yielded poor results. Modes with large coupling to the control system are not 
modeled when natural frequency 1s the method of truncation. Instead modes with very 
little coupling to the contro] svstem are modeled and in attempting to control these 
modes very large control torques are generated causing large excitations in the strongly 
coupled but unmodeled modes as seen in Figures 9 through 14. 

The excitation of each mode depends on two factors: the natural frequency which 
determines the damping and the amplitude of the mode shape which determines how 
much excitation 1s received by each mode. The amplitude of the mode shape is the 
dominate factor which is evident from the open loop responses see Figures 7 and 8. 

Truncating based on the open loop response, 1.e., the Karhunen-Loeve model also 
vields poor results except when it is developed from the open loop response for the case 
of the disturbance torques being due to actuator noise, see Figure 8 on page 25. When 
truncating modes based on the open loop response, the coupling from the noise input 
is the donunate factor in selecting the modes to include in the model. In this case modes 
with large open loop excitation are controlled, but the control coupling still dominates 
yielding large excitation in the unmodeled modes (see Figures 15 through 17). When the 
disturbance torques are due to actuator noise, the modes with the largest open loop 
excitation are also the modes with the largest control coupling and the control system 
works well (see Figures 18 and 19). 

The KL model based on the open loop response for disturbance torques due to 
actuator noise works well when applied to the space station with a disturbance torque 
applied at another location. As seen in Figures 20 through 22, the control system drives 
the modeled modes close to zero without exciting other modes to a significant degree. 
The problem with this configuration is that modes with large open loop excitation are 


not necessarily modeled and may be left excited by the control system. 
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B. DESIGN PROCEDURE 

This design procedure is based on the above conclusions and is recommended for 
controlling vibration in LSS. First, the reduced order control model should be based on 
the KL mode shapes determined from the open loop response for disturbance torques 
due to actuator noise. Next, check the open loop response for disturbance torques at 
locations on the structure where noise inputs are most likely, e.g., the alpha joint. 
Modes strongly coupled to the disturbance input must be included in the KL model. 
This will determine the size of the KL (reduced order) model. If, after all this is done, 
the results are still unacceptable, then adding additional control elements should be 
considered. The placement of the additional control element(s) is determined by the 
node(s) with the largest modal amplitude(s) for the modes that remain excited by the 


initial control system. The process is repeated until an acceptable control is achieved. 
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APPENDIX A. SPACE STRUCTURE SIMULATION PROGRAM 


This program simulates the dual keel space station described in Chapter | by im- 
plementing the model described in Chapter IV. The control described in Chapter IV is 
also simulated. For a detailed explanation of this program and its development see [Ref. 
J}. 


Sevededevevetetevedededededededetedecedtedesetetehi ede tckkiekdietcleckedcdek vevesesesvetekiedcesvedetevck ts 
Wawona Space Structure Simulation Program Tere rese re 
sonenenes < Sceeeeses 


ve ste Je awlee $, Cd ats ale al ale ale 
ey en ee eueeen en 


BY William J. Preston 


ale af ale =‘ ste ste ste steeles ed ‘. se fan o's = aloe = euton'entouts 
PTET Eee Te TEES TOS eae ove DDT ETE TESTES TET TELE RTT TOT TERE TOTP TERETE TET 


Sesteves's 38 se; seve oe we ste ale ale we eves‘: ves'e veves'e aN ee Sevesevoves's ves". we ss * . sieves‘ ¥ sles’ we vey: we sews Seveve fe ute te weve i on 


Cece VaR Brea RAT LORS Wee TT 


sestes'e ale s'e =". evento ~s 3's este s'e af eute a! te wfos) we ses ate s'¢ wes gates ate =‘ ewe at fe ate gteatontenteclonte ate ate afealont ate ale ute afoutertontontontonton'sntoatoutsatentesntonls 
7Cee Ss Cs oe CN EN EN CD EVEN CN ON ER EV CV EC EL CRED EN ENED EN ED EV CLEVER EN EL CL EN OVEN 


GC) 2°61. €) C2) C2 Canta C269" ta 


EXTERNAL EXCMS ,RICDSD 

CHARACTER*6 NAM 

CHARACTER*1 AGAIN,CORECT , RAGAIN 

INTEGER NODE ,MODE ,KQ,EMODE, SMODE 

INTEGER CT,CF,KA, LOOP, PRNT,MODAL,V,COUNT, PRNTG 

INTEGER NF ,NG,NH,NZ,1I,K,M,CrAaADJ 

Noe IPYS( 100) , ITYPE( 200) 

we s'esten' slevlecta a's slo nlenlontentoe te wleals oleate aleate ate steatoetoate a! ate sfecte sie stesle steslecte cle sie 

REAL TOTCST /RMODEN(7: 100) 

ane © cOSW1T, See COST, ele ENERGY , RM 

Saclaster vestesies ‘estewtecten eslevgclesteses io aewtededecied eslsalectente ete ee 

REAL RTOTAL -RMODEN( 7: 100) 

REAL"8 PHII(2,2,100) ,GAMMA( 2,100) ,EGT,GMA,WN,W1,X1T,X2T 
REAL*8 PHI( 188,188) ,B(188,3) ,BN(188,3) ,R(3,3) ,RR(3,3) 
REAL**8 RINV(3,3),RRINV(3,3),X1(7: 100) ,X2(7: 100) ,MODEN( 7: 100) 
REAL*8 COSW1T,SINW1T,COST,CNTCST,ENERGY , TOTCST ,RM 
REAL*8 TCX,TCY,TCZ,DAMP,SAMPT,PI,SUM1,SUM2,SUM3,SUMC 
REAL*8 TNX,INY,INZ, IMPX, IMPY,IMPZ:> THPUSA Sto ee eo 
REAL LAMA( 100) , UGVEX(100,3),RNODE, RMODE ,MIN, TIME, SAMPTM 
REAL UG69(100,3),UG23(100,3),UG55(100, 3) 

REAL*8 H(100,100),G(100,100),L(3,100) ,BT(3, 100) 

REAL*8 Z( 200,200) ,W( 200,200) ,ER(200) ,F( 100, 100) ,EI( 200) 
REAL*8 SCALE( 200) , TEMP( 100,3),TEMP1( 3,100) ,WORK( 100) 
tervervedede dete te dete de Teese te Fete Fete se Ne Ve TNE FEDER FC TCT VF Ve TCT CTT Fe Te Fe Fe Fe TTT ET CTH KKK RRR 


GGG 


Jevetvetevetetetetetevesevededetetedckekictededetededesetekletedevcdedetede deter sek kevesesevetevete dete K tote HK 


detedetete VARIABLE DEFINITIONS fereveeae 
Kictedevededetedsetetedevetedestesdedevevededetesetesevetededetcsetesdescdetkvededetet tevetevetete te kevese te Hei ve 


C2? C2°G2- Gia) 
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So) —) — let etetet@) oR elete | (> (Pe ele elt el elele eee (eS) ei Pele eletelevoteken ene heen en® 


LAMA = VECTOR OF THE SQUARE OF THE NATURAL FREQUENCIES 
UGVEX = MODE POSITONS AND SLOPES OF THE NODAL POINTS 

PHI = STATE TRANSITION MATRICIES FOR EACH MODE 

GAMMA = INPUT TRANSITION MATRIA 

A = DIAGONAL MATRIX CONSISTING OF PHI 

B = INPUT MATRIX OF GAMMA AND CONTROL NODE SLOPES 

BN = NOISE INPUT MATRIX OF GAMMA AND NOISE NODE SLOPES 

DAMP = DAMPING FACTOR 

SAMPT = SAMPLING TIME 

IMPLSE = IMPULSE INPUT FUNCTION 

TCX, TCY, TCZ = CONTROL TORQUE VALUES 

IMPX, IMPY, IMPZ = AXIS IMPULSE NOISE VALUES 

ENERGY = SYSTEM ENERGY COST VALUE FOR A GIVEN POINT IN TIME 
CNTCST = SYSTEM CONTROL COST VALUE FOR A GIVEN POINT IN TIME 
COST = TOTAL SYSTEM COST VALUE FOR A GIVEN POINT IN TIME 
TOTCST = SYSTEM COST SUMMED OVER ALL TIME 

MIN = NUMBER OF MINUTES SYSTEM WILL BE OBSERVED 


Mevevelodesedevetese SAMPLE OF SPACEN EXEC FILE sevetedesetesesesetese sete I 


THIS FILE MUST BEGIN IN COLUMN 1 AND RUN WITH THE FOLLOWING 
SEQUENCE FOR THE INITIAL RUN OF THE PROGRAM: 


FORTVS2 SPACEN (COMPILES PROGRAM) 
SPACEN (LOADS AND RUNS PROGRAM) 


SUBSEQUENT PROGRAM RUNS CAN ELIMINATE ''FORTVS2 SPACEN" IF NO 
CHANGES HAVE BEEN MADE TO THE PROGRAM. 


CP DEF STOR 2M 

FI 4 DISK KLAMA OUTPUT B (PERM 

FI 30 DISK X1 OUTPUT A (RECFM F BLOCK 80 PERM 

FI 31 DISK MODENG SPACEN A (RECFM F BLOCK 80 PERM 
FI 32 DISK TORQUE OUTPUT A (RECFM F BLOCK 80 PERM 
Pi 35° DISK ENERGY OUTPUT A (RECFM F BLOCK 80ePERM 
FI 34 DISK MDECST OUTPUT A (RECFM F BLOCK 80 PERM 
FI 35 DISK COUNT INPUT A (RECFM F BLOCK 80 PERM 
FI 40 DISK UTILITY DATA A (RECFM F BLOCK 80 PER 
FI 41 DISK RUN DATA A (RECFM F BLOCK 80 PERM 

FI 42 DISK KUG69 OUTPUT A (RECFM F BLOCK 80 PERM 
FI 43 DISK KUG23 OUTPUT A (RECFM F BLOCK 80 PERM 
FI 44 DISK KUG55 OUTPUT A (RECFM F BLOCK 80 PERM 
LOAD SPACEN 

START * NOXUFLOW 


Sevetevetededededetevededtedevedededesedededkdesvedcdedetesekededesesdesededesvesededededededcdevevedededededesdesdedesevede 


PI = 4.0D0 * ATAN(1. ODO) 


SAMPT = 0.0 
DAMP = 0.0 
MODAL = 0 
IMPX = 0. ODO 
IMPY = 0. ODO 
IMPZ = 0. 0D0 
NF = 100 
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2 2 Ceo CavG Ga) 


C10) Gi OG) GO) Cy CCG C) 


Gainey I C) Ca) 


Co) E 


100 
100 
NZ = 200 


Severeseslevevesicsk ve ve 


ee a 
4 G) 
tI 


wee NUMBER OF MINUTES THE SYSTEM WILL BE OBSERVED = **%%-% 


MIN = 120.0 


Sesesedededk sevesesevedesde dese ve deveveveseve verve sesesese ve sevevesesesc ve sea sede vesesevese vevesesesesc sesesesesccese sk 
HAWK SET LENGTH OF MODAL MODEL Hewes 


Feveskvestcocslesvlcseveslcve sesesleak vkesdeveske tevesesevevesesesevese severe ve se sesevese ses se seseseseaks cts ae mae ee 


MODAL = 56 


Seveskucorsestestesesteseste tests alte seated steototstestcotse cece asst sk at sisestestaksestseseste st ses ctostese sect seat te se ak oe 


sesesesext READ LAMA MATRIX sevesteccve 
Seovsededesesese ve seveseveseuesededeveskksevesledesevesevevencvescsevedesesedesesesevedcscsdeseseseveveveve ses seve ase 


READ(4,1001) NAM 
READ(4,1002)( LAMA(I) , I=1, 100) 


scoicsesevcvesescsesteslcvevesesesevevevevescsescovsesesesesesesesesesesesecescstesestcsecestesestese score ses arse seslesesess 


Seec sess lt SCREEN INTERACTION Sevevesieve 
sevievesecoveseslecesionesesesesesesevesesevevesedeveveseveseseveve veve vevevese vevesesesevese ses vevevevese te ve veses eres 
skeotesesesledesiedestesleseste STARTING MODE NUMBER Keieseseveslestcsescveslestcslesc 
SMODE = 7 


Pipe lyse hse isviscigeiochaeleelerisetseischs NUMBER OF MODES TO SCAN Sesewesesesevecwes eee 


MODE =5 

EMODE = SMODE + MODE - 1 

Seon yCeCue eC crane NOISE INPUT POSITION Seed sec ec eae eC eC dads sees 
NODE = 55 

AXIS = 1 

rlarverleheriyricr oriveiirlse ipeiseks R MATRIX VALUE HRaevesesesevesesese dese weve ewe 
RM = 1E-12 


AiededeeI O72 FOR FULL SAMPLING TIME °.05 FOR REDUCED ®t 


SAMPT = 0.01 

We Tesese He Tete Te sea ce 6 ae ae DAMPING FACTOR Sivesese tease sessed week 
DAMP = 0.001D00 

DO 75 I = 1,100 


READ( 42,1040) (UG69(1I,K) ,K=1,3) 
READ( 43,1040) (UG23(1,K),K=1,3) 


goon 


Cig ©) 


gqaqgdnsi 


READ( 44,1040) (UG55(1,K) ,K=1,3) 
CONTINUE 
Sesesevedeseseccceseseveseseve BEGIN RUN 


DO 505 M = 1,3 


WRITE (41,700) SMODE 
WRITE (41,701) MODE 
WRITE (41,706) EMODE 
WRITE (41,702) NODE 
WRITE (41,703) RM 
WRITE (41,704) SAMPT 
WRITE (41,705) DAMP 
WRITE (41,707) MIN 
WRITE (41,708) MODAL 


Here Herere vedere NOISE AXIS INPUT AND LOCATION 


IFC AXIS. EQ. 1) THEN 
IMPX = 1.0D0/SAMPT 
ELSEIF( AXIS. EQ. 2) THEN 
IMPY = 1. ODO/SAMPT 
ELSEIF(C AXIS. EQ. 3)THEN 
IMPZ = 1.0D0/SAMPT 
ELSEIF( AXIS. EQ. 4) THEN 


IMPX = 1. ODO/SAMPT 
IMPY = 1. 0D0/SAMPT 
IMPZ = 1.0D0/SAMPT 
ENDIF 
COUNT = QO 
devedevededededededed INITIALIZE MATRICIES 


DO 40 I = 1,188 
DO 45 J = 1,188 
Poa ( 1,3) = 0.0 
CONTINUE 
CONTINUE 


DO 60 I = 1,188 
DO 65 J = 1,3 
Pains) = 0.0 
mer. J) = 0.0 
CONTINUE 
CONTINUE 


DO 70 K=7,1 
paneK i) = 0. 
X2(K) = 0. 
MODEN(K) = 0.0 
RMODEN(K) = 0.0 

CONTINUE 


00 
0 
0 


Jesesescsescoeseseosvevesicve ve 


Seoesesevesesesecese 


sesevesestskeocoksestesess 


fn ata ent, 


sevevese 


Teorcaedeakatacacaedevede sede kak deakokakoatseotcukcskokokevedeskdkcakdeacakcdescvkeakaksesk sk ke kskacacatcsecke seve tee ve ae te ok 


wevedodeds BEGIN MAIN PROGRAM 
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Weslo evo 
ev evar eve 


2 2° @ 


dederesere ESTABLISH A, B AND B'NOISE" MATRICIES dedeseders 


Yeolsvevescdesededededesevededescdetedededkvdevedesededededecdedevedetescdededtedededcdedesededcsdevestesesdeslcdeseslevestestese 


DO 600 I = SMODE,MODAL 
WN = DBLE(SQRT( LAMA(I))) 
GMA = DAMP*WN/2. 0 
EGT = DEXP( -GMA*SAMPT) 
W1 = DSORT( (WN***2)-(GMA**2) ) 
COSW1T = DCOS(W1*SAMPT) 
SINW1IT = DSIN(W1*SAMPT) 


IFCWN. EQ. 0) THEN 


PHII(1,1,I) = EGT*COSW1T 
PHII(1,2,I) = SAMPT 
PHII(2,1,1) = 0 
PHII(2,2,1) = EGT*COSWI1T 
GAMMA(1,I) = 0 
GAMMA(2,I) = 0 
ELSE 
PHII(1,1,1) = EGT*(COSW1T + (GMA*(W1**(-1)))*SINW1T) 
PHII(1,2,1) = (W1**(-1))*EGT*SINW1T 
PHII(2,1,1) = -(WNe*2)*(W1**( -1))*EGT*SINW1T 
PHII(2,2,1) = EGT*(COSW1T - (GMA*(W1**(-1)))*SINWIT) 
GAMMA(1,1) = (WN%¥***(-2))**(1. ODO-EGT*( COSW1T+(GMA/W1)*SINW1IT)) 
GAMMA(2,I) = (W1%*(-1))*EGT*SINW1T 
ENDIF 
CONTINUE 
ve 


DO 610 K = SMODE,MODAL 


PHI(V,V) = PHII(1,1,XK) 
PHI(V,V+1) = PHII(1,2,K) 
PHI(V+1,V) = PHII(2,1,K) 
PHI(V+1,V+1) = PHII(2,2,K) 


B(V,1) = GAMMA(1,K)*DBLE(UG69(K,1)) 
B(V,2) = GAMMA(1,K)*DBLE(UG69(K, 2)) 
B(V,3) = GAMMA(1,K)*DBLE(UG69(K, 3)) 


B(V+1,1) = GAMMA(2,K)*DBLE(UG69(K,1)) 
BC V+1,2) = GAMMA(2,K)*DBLE(UG69(K,2)) 
B(V+1,3) = GAMMA(2,K)*DBLE(UG69(K,3)) 
V = V+2 

CONTINUE 


DO 605 I = 1,100 
UGVEX(I,1) = 0.0 
UGVEACI,2) = 0.0 
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605 


C 


Kho 
= 


MQOoQUNI ONAN 
‘Tal 
© 


55 


66 
61 


UGVEX( 1,3) = 0.0 
IFCNODE. EQ. 23) THEN 


UGVEX(I,1) = UG23(I,1) 
UGVEX(I,2) = UG23(I,2) 
UGVEX(1I,3) = UG23(1,3) 

ELSEIF(NODE. EQ. 55) THEN 
UGVEX(I,1) = UG55(I,1) 
UGVEX(I,2) = UG55(I,2) 
UGVEX(I,3) = UG55(I,3) 

ELSEIF(NODE. EQ. 69) THEN 
UGVEX(I,1) = UG69(I,1) 
UGVEX(1I,2) = UG69(I,2) 
UGVEX(I,3) = UG69(I,3) 

ENDIF 

CONTINUE 


V=1 
DO 620 K = SMODE ,MODAL 


BNCV,1) = GAMMAC1,K)*DBLECUGVEX(K,1)) 
BNCV,2) = GAMMAC1,K)*DBLE( UGVEX(CK, 2) ) 
BN(CV,3) = GAMMAC1,K)*DBLECUGVEX(K,3)) 
BNC V+1,1) = GAMMA(2,K)*DBLE( UGVEX(K,1) ) 
BNC V+1,2) = GAMMA(2,K)*DBLECUGVEX(K, 2) ) 
BN(V+1,3) = GAMMA( 2 ,K)*DBLE( UGVEX(K,3)) 
V = V+2 
CONTINUE 
CONTINUE 
BeemeCsanaecetaecc ec ae ESTABLISH HS; F AND R MATRICIES Becertac we TCaC wae re 


DO 50 I = 1,NH 
DO 55 J = 1,N 
Hers) 
met J) 
Git) 
TF(I. LE. 3 
fer,J) = 
ENDIF 
CONTINUE 
CONTINUE 


wou i 
wooo; nt 


0 
0 
0 
THEN 
0.0 


Horel I = 1,3 
DO 66 J = 1,3 
R(I,J) = 0.0 
RINV(I,J) = 0.0 
CONTINUE 
CONTINUE 
KQ = 1 
DO 80 K = SMODE,EMODE 
H(KQ,KQ) = DBLE(LAMA(K)) 
H(KO+1,KQt+1) = 1. 0DO 
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© 


mamma 


104 
EOS 


KQ = KQt+2 
CONTINUE 


io=n0 
DO 165) kK eee 

R(K,K) = RM 
RINV(K,K) = 1. 0D0/RM 
CONTINUE 


3 


DO 88 I = 1,2*MODE 
DO 89 J = 1,2*MODE 
F(I,J) = PHI(I,J) 
CONTINUE 
CONTINUE 


wvstevere = «COMPUTE G MATRIX AS NEEDED BY RICDSD SUBR _—*vevevexe 
CALL MATRAN(B, 188, 2**MODE ,3,BT,3) 
CALL MATMUL(B,188, 2**MODE,3,RINV,3,3,TEMP,NH) 


CALL MATMUL( TEMP ,NH, 2“=MODER3 . BI) 3, 2-MODE GanG 


wlaalasuls u's waafon'satsaisuisuloniontsetcu'ostasn'sclienfsnioniantian'anou'eeiesieniassionsuonieniouanseenlaniontontsniouants We nie wlaales wontentoneniacteutentiousulon'ou'e 
Peed ER CR RES CR CR Ed FV EDM CR FR EX FL FT ERED ELF do FL FX FR EL EDS GO FL CR FL Ge FE FE FEEL IR GL CO Fo #% #% 38 @¥ FX EV ER FL Fo Ro FR FR EN GO. FS CoE ER FL #% 


BEGIN RICCATI GAIN CALCULATIONS 


wleates wJwaleesalsasioniastoniontonionionioutonionanisuenisn ou'snouaevionieniononssisuenianionioslieneuansuenianiontenianosionauliouonioucuianian! wlontes wlontontiou's 
ER FR EU ER EL ER FR TR EH ER EL EM ED FY ER FT ELEN EL FRED FR EY FY ER CO EL ER PE GE FL EL ER FER EL GO FH FL ER FR ER ER FR CXF EL eS ED HX ER FL FD an OV GX ed. FX FX FD 


CALL RICDSD(NF,NG,NH,NZ,2**MODE, 4**MODE,F,G,H,Z,W,ER,EI,WORK, 
SCALE, ITYPE, IPVS) 


WRITE(41,*) ' ' 
WRITE (41,130) 
FORMAT (/' THE CLOSED LOOP EIGENVALUES ARE: '/) 
DO 140 I = 1,2**MODE 
WRITE (41,*) ER(1I),EI(I) 
CONTINUE 
WRITE (41,150) WORK(1) 
FORMAT (/' CONDITION ESTIMATE IS: ' ,D26. 18) 


Sen aia mi sie ac aie ce he ee COMPUTE GAIN MATRIX L SOF yey omen en eae secs sce e 
CALL MATMUL(BT,3,3,2%*MODE,H,NH, 2**MODE, TEMP1, 3) 
CALL MATMUL(TEMP1,3,3,2*MODE,B,188,3,RR,3) 
DO 103 I= 1,3 
DO 104 J=1,3 
RR(I,J)=R(I,J)+RR(I,J) 
CONTINUE 
CONTINUE 
CALL DLINDS(3,RR,3,RRINV,3) 


CALL MATMUL(RRINV,3,3,3,TEMP1,3,2*MODE , BT, 3) 
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CrArieMAIMUEGRT,3,3,2*MODE,PHI ,188,2*MODE,L,3) 


C 
WRITE(41,*) ' ' 
WRITE(41,*) ' GAIN MATRIX L '' 
WRITE(41,*) ° ' 
WRITE(41,*) ' ROW 1 ROW 2 ROW 3' 


DO 9155 I=1,2*MODE 
WRIBE G4. 1040 )m@ECJ.1) ,J=1,3) 
9155 CONTINUE 
WRITE(41,*) ' ' 


C 
CG Jetedededed tessvededede dedesvevededevedede te teve de Fede te FoF Te Te EK Fe eFC TEN TCT Fe He Ne Fe Fete Tete se Kee 
c dederedere COMPUTATION OF TORQUES AND COSTS derevesede 
Cc Sc aanuaens amar 1622s 26 ay Tombs 26 20 90 26 Je cna C eC an ae aC 2s aC eC IE Oe OC TC 76 DC OE Oe SE De OC DEE IE Oe TE TE ae IDC Ie TOC ee ae Ie Oe 
C 
9000 COUNT = 0 
TOTCST = 0. ODO 
TIME = 0.0 
C 
G Bere re seve Te ae se 76 BCI Fee F076 Ie I0 IC IT IE IE Bh IE IE EI NE IEE DET TCDS IE IEDC EME IE ETE IE IE TE IED IE TE TENE MEME DE WIE DE DE DE DEVE IE WE VE 
C vekeveere SETS LOOP FOR THE NUMBER OF ITERATIONS NECESSARY — *vtveves 
c yee ~=TO OBSERVE THE SYSTEM FOR DESIRED LENGTH OF TIME ‘reves 
Cc Serkedededevededevevevedevedededevevesedevededeve Tevevedede deveve Seve k ve Fete te Seve deve ve He sete de fe te se sene Joveve sede ves eve ve 
C 
LOOP = INTC CMIN*60. 0)/SAMPT) 
PRNT = INTC CCMIN*60. 0) /SAMPT) /300. 0) 
PRNTG = INTC CCMIN*60. 0)/SAMPT)/1000. 0) 
C 
DO 200 N = Q, LOOP 
TIME = DBLE(N)**SAMPT 
C 
IF(N. EQ. 0)THEN 
IMPLSX = IMPX 
IMPLSY = IMPY 
IMPLSZ = IMPZ 
ELSE 
IMPLSX = QO. ODO 
IMPLSY = Q. ODO 
IMPLSZ = Q. ODO 
ENDIF 
C 
Cc >No cos ni th arb tdi ee Gi bie ie Siok Or ie Sica Sk So OCICS CSIC ICI aici ba Oi oi ok Gini Giro Oi ieGaei inbird bir id 
C devededesevededsdeds CONTROL TORQUE EQUATIONS devededededodededsseds 
GC POE oo eee ee ee tock cece eae cS ca Oe cece ee IN Se eae Ee RSE eT TE RUC TR SATE TCT OTTO TETRIS ESE 
C 
SUM1 = QO. ODO 
SUM2 = Q. ODO 
SUM3 = 0. ODO 
C 
DO 210 CT = 1, MODE 
CTADJ = CT + (SMODE - 1) 
SUM1 = SUM1 + L(1,2*CT-1)*X1(CTADJ) + L(1,2*CT)*X2(CTADJ) 
SUM2 = SUM2 + L(2,2*CT~-1)*X1(CTADJ) + L(2,2*CT)*X2(CTADJ) 
SUM3 = SUM3 + L(3,2*CT-1)*X1(CTADJ) + L(3,2*CT)*X2(CTADJ) 
210 CONTINUE 
C 
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Pe Caeterer re C7 Ceca o) C2) OO) MOI OO Ome OI Oe Om OE CoC CG) "C0 GC) Ca 


aC). G2) C2 


Fervdereters deve se He Hye EVALUATION AT ZERO CONTROL Tedededevededededese deve 


TECH. SQ. Gnas 


TCX = Q. ODO 
TCY = 0. ODO 
TCZ = QO. ODO 
ELSE 
TCX = SUMI*( -13050) 
TCY = SUM2*(-1. ODO) 
TCZ = SUM3*(-1. ODO) 
ENDIF 


IF(N. EQ. 0) THEN 
WRITE (32,*) ‘IMPULSE X AXIS, IMPULSE Y AXIS, IMPULSE Z AXIS' 
WRIMGE (.324%)e oe 
WRITE (32,1040) IMPLSX, IMPLSY, IMPLSZ 
WRITE SC32 °° ae ae 
WRITE (32,%) 'CONTR@E. TORQUBS TCX, TCY, TGZ' 
WRITE 60325" a 
ENDIF 


IF (N. LE.50) THEN 
WRITE (32,2000) TIMES TCx, Tey, 1CzZ 
ENDIF 


IF (MOD(N,PRNTG).EQ. 0) THEN 
WRITE (32,2000) Tithe [Cre ren, ez 
ENDIF 


IF (N. LE. 20) THEN 
WRITE( 30,1035) TIME,X1(1),X1(5),X1(10),X1(15),X1(20) 
WRITE( 31,1035) TIME,X2(1),X2(5) ,X2( 100), X2@s eee 
ENDIF 
IF (M.EQ. 1) THEN 
IF (MOD(N,PRNTG).EQ.0) THEN 
WRITE( 30,1036) TIME,X1(9),X1(10),X1(11),X1(12),X1(13) ,X1(14) 
WRITE( 30,1036) TIME,X1(1),X1(4) ,X1(14) ,X1( 24) ,X1(34) ,X1(44) 
WRITE( 31,1035) TIME,X2(1),X2(5),X2(10),X2(15),X2(20) 
ENDIF 
ELSE 
IF (MOD(N,PRNT).EQ.0) THEN 
WRITE( 30,1036) TIME,X1(1),X1(4) ,X1(24),X1(44) ,X1(74) ,X1(94) 
WRITE( 30,1036) TIME,X1(1),X1(4) ,X1(14),X1(24) ,X1(34) ,X1(44) 
WRITE(31,1035) TIME,X2(1),X2(5),X2(10),X2(15) ,X2(20) 
ENDIF 
ENDIF 


IF (MGDGNSERN@G).EO. 0) THEN 
COUNT = COUNT+1 
ENDIF 


Krcsetevestedevetevedetescakeseseukesescsestcoescscskcskcslesestcokesescsicsestesicaicakcskcseskestestcse ak siesta se ae ales se ak ve ae se ve ok 


Weveredederedede SYSTEM COST FUNCTION CALCULATION Hee 
Ticdevevedededededtededetevedevedesesetevevesedeckestcdcskedevevesedeseakeseveveve teak vested sete vkedeakakeskcoke devesese ve seat 
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Gd 
>) 


C202 C7 Corre? C202 C2 


C2)C2 C2 Sa a 


235 


666 


SUMC = QO. ODO 
ENERGY = 0. ODO 
CNTCST = Q. ODO 
COST = 0. ODO 


DO 230 CF = 7,MODAL 
MODEN(CF) = MODEN(CF)+(X1(CF)**2)*LAMA( CF)+X2(CF)**2 
SUMC = SUMC+(X1(CF)*2)*LAMA( CF )+X2(CF)**2 

CONTINUE 


ENERGY = SUMC 

CNICST = (TCX**2)*RM+( TCY** 2 )*RM+( TC Z** 2 ) *RM 
COST = ENERGY + CNTCST 

TOTCST = TOTCST + COST 


IF (MOD(N,PRNTG).EQ.0) THEN 
feiss 2000) iin iNERGY ,CNICST,COST 
ENDIF 
PRON. GE. (EROOPR-50)) THEN 
Wiki 34, 5002) MODE, TOICST 
ENDIE 


Sescsicdeacasesescsvscsestcsedesestckseseslesesleskslevesiestcsvsleskskevleseseskseslecsesesevescccsesesleskcovveseseslestesleseste sk se 
Fovedeaede detec STATE UPDATE EQUATIONS deaededkdodededesede 


aa 


SES Se CIS Fe OE CTR I IE TE Ie CSC OS See re ae oe sete te sesesesestseveseseskeseck dese seseskeseskvesesesesevess sess 


DO 220 KA = 7,MODAL 


K = KA-6 
X1T=PHII(1,1,KA)*X1(KA)+PHII(1,2,KA)*X2(KA)+B((2*K-1) ,1)*#TCX+ 
+ B((2*K-1) ,2)*TCY+B((2**K-1) , 3)*TCZ+BN((2*K-1),1)*IMPLSX+ 
BN( (2*K-1) ,2)*IMPLSY+BN((2**K-1),3)*IMPLSZ 


X2T=PHII(2,1,KA)*X1(KA)+PHII(2,2,KA)*X2(KA)+B(C2"K, 1)*TCX+ 
B( 2*K, 2)*TCY+B( 2*K , 3)*TCZ+BN( 2*K, 1)*IMPLSX+BN( 2*K ,2)** 


+ IMPLSY+BN( 2*K,3)*IMPLSZ 


X1(KA) 
X2(KA) 


on 
X2T 


CONTINUE 
CONTINUE 


WRITE (35,3000) COUNT 
RomAL = TORCST 

WRITE (34,3002) MODE,RTOTAL 
WRITE (34,3002) MODE,TOTCST 


DO 235 K = 7,MODAL 

RMODEN(K) = MODEN(K) 

WRITE (31,3002) K, RMODEN(K) 
CONTINUE 


CONTINUE 


mM 


2001 
3000 


IP @iREO SiN 


MODE=10 
EMODE=16 
ELSEIF(M. EQ. 2)THEN 
MODE=15 
EMODE=21 
ELSEIF(M. EQ. 2)THEN 
MODE=20 
EMODE=26 
ELSEIF(M. EQ. 2) THEN 
MODE=25 
EMODE=31 
ELSEIF(M. EQ. 3)THEN 
MODE=30 
EMODE=36 
ELSEIF(M. EQ. 4)THEN 
MODE=35 
EMODE=41 
IF(M. EQ. 1)THEN 
MODE=40 
EMODE=46 
ELSEIF(M. EQ. 2) THEN 
MODE=45 
EMODE=51 
ELSEIF(M. EQ. 3)THEN 
MODE=50 
EMODE=56 
ENDIF 
CONTINUE 
tesededesdeskeslescdesedlese dk okvacakcslesedevledescsdkeslesesccdededcsdedescacsdesedededesd dedevledeacdcakcsesccedeakeske sede ak dese ve se 
Veseveces'e FORMAT STATEMENTS Sevesloveste 
aedesevedededenededededeseseteskeakdededesevedededededesededevededededcdcdededevededededesdcdescdededesdededesdescdededese 
FORMAT (' ','STARTING MODE NUMBER: ‘',12) 
FORMAT (' ',' NUMBER OF MODES SCANNED: ‘' ,12) 
FORMAT (' ','NOISE INPUT NODE: ',I3) 
FORMAT (' ', ‘INITIAL R VALUE: ‘',E12.4) 
FORMAT (' ','SAMPLING TIME: ‘',E12.4) 
FORMAT (' ','DAMPING FACTOR: ‘',E12. 4) 
FORMAT (' ','LAST CONTROLLED MODE: ‘',I2) 
FORMAT (' ', ‘OBSERVATION TIME: ',F5.1,' MINUTES’) 
FORMAT (' ','SIZE OF MODAL MODEL: ',13,' MODES’) 
FORMAT( 1X, A6) 
FORMAT( 1X,8E15. 8) 
FORMAT(1X,//) 
FORMAT(1X,60X,E11. 5) 
FORMAT(1X,////) 
FORMAT(A1) 


FORMAT(' ',F7.2,2X,5(E12. 6,2X)) 
FORMAIC” " ,F7.2,1X%.6¢2a175 lop 
FORMAT(' ' ,3(E15.8,5X)) 
FORMAT(' ',4(E12.6,2X)) 

FORMAT( 1X,F7. 2,3X,3(E15. 8,3X)) 
FORMAT ( '@' )BSepkMiseesa) 

FORMAT( 14) 
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3001 
3002 


299 


Care) Ca C2: ©) 


670 


660 
650 


C21@ Cr C2 


690 
680 


FORMAT(F7. 2,2X,E12. 5) 
FORMAT(I3,2X,E12. 5) 


ShOly 
END 


sesesteclesicskoesesesesesevescdesevetedesoveveveseveslesese se desesesesevededevevenesedevesic ak sleslesic ves sieve seskskeskeske desk 


SUBROUTINE TO MATRIX MULTIPLY 


seveveavcesescatcevesicsesestestecslesesesdesesedeskscvlcskcsvseakesescvedovesenesentsicvleskcslest akvsededeskcslevesesevestcskesecese ve 


SUEROULINE MATMUL(UML, LDL, RI,Ci M2,L02,C2,MP,LD3) 
TNOEGERSR1,Ci.62, LD, LP2YLEDS 
fe Mel LDL, 1) ZC? ey MPC LD3 ,1) ,SUM 


DO 650 I = 1,R1 
DO 660 J = 1,C2 
SUM = 0. ODO 
DO 670 K = 1,C1 
SUM = SUM+M1(1I,K)*M2(K,J) 


CONTINUE 
MP(I,J) = SUM 
CONTINUE 
CONTINUE 
RETURN 
END 


Sesdedevkededestesiovevesoveksescdodksksdlevcdesesedcveskultevevededesleslesteslesievevlecsedecdicsesienk sleskevicsiesestesleskesicdestese ak sts'c 


SUBROUTINE TO TRANSPOSE A MATRIX 


sesicsiesiesicsesleseslesesiesesesesedeslesievesieseslesesesesesedevevesesedeseslesleneslestesicvkcslescvestedeseseveveslescvieskedeslesieste sie 


SUBROUTINE MATRANCMX,LDX,R1,C1,MT, LDT) 
PIPEGERSRI,C1,1,J, DX, EDT 
REAL*8 MX(LDX,1),MTCLDT,1) 


DO 680 I = 1,R1 
DO 690 J = 1,C1 
MT(J,I) = MX(I,J) 
CONTINUE 
CONTINUE 
RETURN 
END 


APPENDIX B. KARHUNEN-LOEVE MODE PROGRAM 


This program computes the Karhunen-Loeve mode shapes as discussed in Chapter 
IV. Subroutine GENRAA 1s called to generate the covariance matrix of the modal co- 
ordinates, RAA=E[a(t)at(t)] . One subroutine from the IMSL math library 
(DMRRRR) 1s used for matrix multiplication and two EISPACK subroutines are used 
to compute the eigenvalues and eigenvectors of the covariance matrix. 

The subroutine GENRAA as noted above generates the covariance matrix, RAA, 
of the modal coordinates of the space station discussed in Chapter I. The algorithm: 


¢ Computes the natural frequency, damping coefficient, y, and damped frequency of 
each mode 


e Computes the variances and covariances based on these parameters as given by 
equations 19, 20, 21 


The KL mode shapes are found by solving the eigenvalue problem, equation 12. 


RAA*Q¢ =1¢6 (35) 


where the matrix Q 1s: 


Q= ing [ | i 7 - i 7 (3G) 
and 
e / 1s the eigenvalue of the product, RAA * Q 
e represents the orthonormal eigenvectors 
RAA * Q is not a symmetric matrix so the eigenvalue problem is changed to: 
[Q'* RAA Q?] [Q'* 6] =4[Q"" 6] (37) 
iret 
n=Q'*pand@=Q'’n (38) 
and finally: 
[Q'* RAA QU] n=An (39) 
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The EISPACK routines for symmetric matrices are used to solve the eigenvalue prob- 
lem, equation 39, and equation 38 1s used to compute the KL mode shapes. 
This program is run by executing the following EXEC file: 

¢ FORTVS2 KLMODES2 

= ies ISK THESIS INPUT BPE 

e FI6 DISK RAA OUTPUT T (PERM 

* Els DISK EIGVEC OUTPUT (PERM 

rl we DISK EIGVAL OUTPUT T (PERM 

e FI] 09 DISK EIGEN ERROR B (PERM 

e FI 12 DISK NATMODES DOO! B (PERM 

e FI 13 DISK KLAMAQ OUTPLT B (PERM 

e FI 14 DISK RAA OUTPUT (PERM 

e FI 17 DISK NERROO] DATA A (PERM 

ar DEareStOR 2M 

© EXEC MATHPACK EISPACK 

meEXECTDISK 5 DIS 

® LOAD KLMODES? (START 


ee KLMODE 
IdekkkkdvdedededekseiekKei Variable Gefin it LONS KKK TS FIT TET eI HIER TE TERT He 


¥ lama = vector of natural frequencies 


ute 
¢% 


< ugvex= matrix of the modal amplitudes and modal slopes * 

* natmod=vector containing the natural mode shapes that correspond ~* 
* to the KL mode shapes a 

* klama=vector of natural frequencies ordered according to th KL 

* mode shapes ¥ 

¥ klmod=vector containing the largest eigenvectors in column of the * 
* ac ieee matrix eg * 


Oe et Re) a i on ee! ao ee) a te en a es) ~ <a awe mia wow ons nlem'o eles wl ef els ol ol ee 
Jeve seve ¥s ve sees Teves sevesieves's se re ss ee Yevs — se se =. ate st Yeves' Jevle ws ves eve eve ses‘ feve eves gues festevs eves sie oe Seve rh ws ws ws ne zu ew a a ae ex ¢d @e ed ws 


CHARACTER*6 NAM 

DOUBLE PRECISION LAMA(100,1),UGVEX(684, 100) ,D(100) ,E(100) 
DOUBLE PRECISION KLMOD(100),Y,N1,Q( 100,100) 

DOUBLE PRECISION GAMA,DAMP,W(50),L(50),M(50),B2(1,50) ,Q2( 100, 100) 
DOUBLE PRECISION K(1,50),A,B,C,COEF,RAA(100, 100) ,APB,APC,TEMP1 
DOUBLE PRECISION TEMP2 ,TEMP3,TEMP4, TEMP5 , TEMP6 ,Q2RAA( 100,100) 
DOUBLE PRECISION EIGVEC( 100,100) ,Z(100, 100) ,RAAQ( 100, 100) 
DOUBLE PRECISION Q2RQ2(100, 100) ,Q2INV(100, 100) ,KLAMA(50) 
INTEGER KK, LL,ZZ,MODR( 100) ,MODC( 100) ,MM,NATMOD( 100) ,NUM 

DO 160 I=1,100 

DO 170 J=1,100 

Q2(1,J)=0. ODOO 

Q2INV(1,J)=0. ODOO 

RAA(I,J)=0. ODOO 


D9 


170 
160 


Q(I,J)=0. ODOO 
CONTINUE 

CONTINUE 

READ(4,1001) NAM 

READ(4,1002) (LAMA(I,1),I=1,100) 
READ(4,1001) NAM 

READ(4,1002) ((UGVEX(I,J),I=1,684) ,J=1,100) 


C NODE IS THE DISTURBANCE INPUT LOCATION  , 
C AND DILOC IS THE COMPUTED LOCATION OF THE MODAL SLOPE IN UGVEX 


NODE=69 
DISLOC = 4+6*(NODE-1) 


C NUMODE IS THE NUMBER OF THE LAST MODE CONSIDERED 


150 


15 


NUMODE = 56 

CALL GENRAA(NUMODE, DISLOC, LAMA ,UGVEX,W,L,M,B2,K,RAA) 
WRITE(16,1005)WC(1),LC1),MC1) , Bo Gu re a 

a, 

tei 

DO 150 J=7 ,NUMODE 

Q2(KK,KK)=DSQRT( LAMA(J,1)) 

U2CLL min — i O00 

Q2INV( KK , KK)=1/Q2( KK, KK) 

G2INVCLL, LL)—l. 0beo 

KK=KK+2 

LL=LL+2 

CONTINUE 

CALL DMRRRR( 100, 100, Q2), 1009100, 100, 92, LOOMMGO, 10050, 16G8 
CALL DMRRRR(100,100,RAA,100,100,100,Q,100,100,100,RAAQ,100) 
DO 15 I=1,100 

WRITE( 14,1007) RAA(I,I) 

WRITE( 12991007) Q2(1,1) 

CONTINUE 

CALL DMRRRR(100,100,Q2,100,100,100,RAA,100,100,100,Q2RAA,100) 
CALL DMRRRR(100,100,Q2RAA,100,100,100,Q2,100,100,100,Q2RQ2,100) 
CALL TRED2(100,100,Q2RQ2,D,E,Z) 

CALL TQL2(100,100,0,E.2, lene 

IF( IERR. NE. 0)THEN 

WRITE(9,*) 'ERROR=' , 16 

PAUSE 

ENDIF 

CALL DMRRRR(100,100,Q2INV,100,100,100,Z,100,100,100,EIGVEC, 100) 
WRITEC6 , 1005), GCRAAGI J) oJ —1 Oe ele) 
WRITE(8,1006)((EIGVEC(I,J) ,J=1,100),I=1,100) 
WRITE(10,1007)(D(1I),I=1,100) 

MM=1 

DO 10 I=100,1,=1 

KLMOD( MM)=0. ODOO 

DON20) J—14 Go 

IF( DABS(KLMOD( MM) ). LT. DABS(EIGVEC(J,1I)))THEN 

KLMOD( MM)=EIGVEC(J,I) 

MODR( MM)=I 

MODC( MM )=J 

Y=(DBLE( MODC(MM)))/2. ODOO 

T1=INT(Y) 

N1=Y-( DBLE(I1)) 

IF( DABS(N1). NE. 0. ODO00)THEN 

NATMOD( MM )=(MODC(MM)4+1)/2 + 6 


56 


20 


10 


30 


180 


ak 
100 


1001 
1002 
1005 
1006 
1007 


100 


iLO 


ELSE 
NATMOD(MM)=MODC(MM)/2 + 6 
ENDIF 
ENDIF 
CONTINUE 
MM=MM+1 
CONTINUE 
DO 30 I=1,100 
WRITE( 12,100) MODR(I),MODC(1),NATMOD(1),KLMOD( I) 
CONTINUE 
Z2=50 
homeo 1=1,100,2 
KLAMA( ZZ )=LAMA( NATMOD( I), 1) 
ZZ=22-1 
CONTINUE 
WRITE(13,1005)(KLAMA(1),I=1,50) 
NUM=7 
DO 11 I=1,99,2 
IF(DABS(KLMOD( I) ). GT. DABS( KLMOD( I+1) ) ) THEN 
WRITE( 17,200) NUM,DSQRT( 2*( 1-DABS(KLMOD(1I)))) 
ELSE 
WRITE(17,200) NUM,DSQRT( 2*( 1-DABS(KLMOD( I+1)))) 
ENDIF 
NUM=NUM+1 
CONTINUE 
FORMAT( 1X, 13,2X,13,2X,13,2X,E15. 8) 
FORMAT( 1X,12,2X,E15. 8) 
FORMAT( 1X, A6) 
FORMAT( 1X, 8E15. 8) 
FORMAT(1X,5E15. 8//) 
FORMAT( 1X,5E15. 8//) 
FORMAT( 1X,5E15. 8//) 
STOP 
END 
SUBROUTINE GENRAA (NUMODE,DISLOC,LAMA,UGVEX,W,L,M,B2,K,RAA) 
DOUBLE PRECISION LAMA(100,1),UGVEX(684, 100) 
DOUBLE PRECISION GAMA, DAMP,W(50),1(50),M(50),B2(1,50) 
DOUBLE PRECISION K(1,50),A4,B,C,COEF,RAA(100,100),APB,APC,TEMP1 
DOUBLE PRECISION TEMP2,TEMP3,TEMP4,TEMP5 , TEMP6 
DAMP=0. 001D00 
GAMA=DAMP/2. ODOO 
DO 100 MODE=7 ,NUMODE 
W(MODE-6 )=DSQRT( LAMA (MODE,1)) 
L(MODE-6 )=W( MODE -6 )**DSQRT( 1. ODO00-GAMA*GAMA) 
M( MODE -6 )=GAMA**W( MODE -6) 
CONTINUE 
DO 110 J=7,NUMODE 
B2(1,J-6)=UGVEX( DISLOC+4**( J-1),J) 
K(1,J-6)=B2(1,J-6)/L(J-6) 
K(1,J-6)=UGVEX(DISLOC,J)/L(J-6) 
CONTINUE 
DO 130 I=1,50 
DO 140 J=1,50 
=M(1)+M(J) 
B=L(I)-L(J) 
C=L(1)+L(J) 


37, 


140 
130 


COEF=K(1,1I)*K(1,J) 

APB=2*( A**A+B**B ) 

APC=2*( A*A+C*C ) 

RAA( 2*I-1,2*J-1)=COEF*(A/APB-A/APC) 
TEMP1=L(J)*B/APB+L(J)*C/APC 

TEMP2=M(J)*A/APB-A*M( J) /APC 

RAAC2* 1-1, 2*J)=COEP*( TEMP Se 2) 

TEMP3=( LC 1I)*LCJ)*A-LCI)*MCJ)*B-LCI)*MC1I)*B+M(1)*MCJ)*A) /APB 
TEMP4=(LC1)*L(J)*A-LCI)*MCJ)*C-LO1T)*MC1)*C-MCT)*MCJ)*A) /APC 
RAA( 2*JT , 2*J )=COEF*( TEMP3+TEMP4 ) 

TEMP5=L(1)*C/APC-L(1)*B/APB 

TEMP6=M(1)*A/APB-MCI)*A/APC 

RAA( 2*I ,2*J-1)=COEF*( TEMP5 -TEMP6) 

CONTINUE 

CONTINUE 

RETURN 

END 
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